Effect of dynamical traps on chaotic transport in a 

meandering jet flow 

M.Yu. Uleysky, M.V. Budyansky, and S.V. Prants 
December 21, 2011 



Abstract 

We continue our study of chaotic mixing 
and transport of passive particles in a simple 
model of a meandering jet flow [Prants, et al, 
Chaos 16, 033117 (2006)]. In the present pa- 
per we study and explain phenomenologically 
a connection between dynamical, topological, 
and statistical properties of chaotic mixing 
and transport in the model flow in terms of 
dynamical traps, singular zones in the phase 
space where particles may spend arbitrary 
long but finite time [Zaslavsky, Phys. D 168— 
169, 292 (2002)]. The transport of passive 
particles is described in terms of lengths and 
durations of zonal flights which are events be- 
tween two successive changes of sign of zonal 
velocity. Some peculiarities of the respec- 
tive probability density functions for short 
flights are proven to be caused by the so- 
called rotational-islands traps connected with 
the boundaries of resonant islands (includ- 
ing the vortex cores) filled with the parti- 
cles moving in the same frame and the saddle 
traps connected with periodic saddle trajec- 
tories. Whereas, the statistics of long flights 



can be explained by the influence of the so- 
called ballistic-islands traps filled with the 
particles moving from a frame to frame. 



1 Introduction 

Chaotic advection of water masses with 
their physical and biochemical charac- 
teristics in quasi two-dimensional geo- 
physical flows in the ocean and atmo- 
sphere can be studied within the frame- 
work of Hamiltonian dynamics. In a re- 
cent paper [1] we have studied chaotic 
mixing and transport of passive parti- 
cles in a simple kinematic model of a 
meandering jet flow motivated by the 
problem of lateral mixing in the west- 
ern boundary currents in the ocean. 
We found all the possible bifurcations 
of advection equations, described the 
structure of the phase space (which is 
the physical space for advected par- 
ticles), and computed some statisti- 
cal characteristics of chaotic transport. 
In the present paper we establish a 
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phenomenological connection between 
dynamical, topological, and statisti- 
cal properties of chaotic transport and 
mixing in the same flow. Speciflc sin- 
gular zones in the phase space where 
particles may spend arbitrary long but 
flnite time (dynamical traps in termi- 
nology by Zaslavsky [2J), are respon- 
sible for anomalous statistical proper- 
ties. The dynamical traps, connected 
with rotational islands and saddle tra- 
jectories, are responsible, mainly, for 
anomalous mixing, whereas those ones, 
connected with ballistic islands — for 
anomalous transport. These dynamical 
traps may have strong impact on trans- 
port and mixing in real geophysical jet 
flows. 



Methods of the theory of dynamical sys- 
tems are actively used to describe advection 
of water (air) masses and their properties in 
the ocean and atmosphere jT] [21 [T31 [151 [ISl 
[26l[3lil|5l[6lEl[I71[T8lC9l[20l[2ll|32]. The 
geophysical jet currents, like the Gulf Stream, 
the Kuroshio, the Antarctic circumpolar cur- 
rent, and others in the ocean and the polar 
night Antarctic jet in the atmosphere, are ro- 
bust structures whose form typically changes 
in space and time in a meander-like way. If 
advected particles rapidly adjust their own 
velocity to that of background flow and do 
not affect the flow properties, then the par- 
ticles are called passive (scalars, tracers, or 
Lagrangian particles) and their equation of 
motion is very simple 

-JT = ^{r,t) (1) 



where r = (x, y, z) and v = (u, f , w) are the 
position vector and the particle velocity at 
a point {x,y,z). If the corresponding Eu- 
lerian velocity field is supposed to be regu- 
lar, the vector Eq. ([1]) in nontrivial cases is a 
set of three nonlinear deterministic differen- 
tial equations whose phase space is a physical 
space of advected particles. It is well known 
from dynamical systems theory that solutions 
of this kind of equations can be chaotic in 
the sense of exponential sensitivity to small 
variations in initial conditions and/or con- 
trol parameters. As to advection equations, 
it was Arnold |i8j who firstly suggested chaos 
in the field lines (and, therefore, in trajecto- 
ries) for a special class of three-dimensional 
stationary flows (so-called ABC flows), and 
this suggestion has been confirmed numer- 
ically by Henon [9]. In the approximation 
of incompressible planar flows, the velocity 
components can be expressed in terms of 
a stream function [TD]: u = —d'^/dy and 
V = d'^/dx. The equations of motion ([T]) in 
a two-dimensional incompressible flow have 
now a Hamiltonian form with the streamfunc- 
tion \& playing the role of a Hamiltonian and 
the coordinates (x, y) of a particle playing the 
role of canonically conjugated variables. 

Advection of passive particles has been 
shown to be chaotic in a number of theoreti- 
cal |a[5l[g[71[Il[Tl[T3l[T5l[l6l[I71|26l[32] 
and laboratory [IHl [ISl [201 [21] models of 
geophysical jet currents. In our recent pa- 
per [Ij (Paper I) we have studied mixing, 
transport, and chaotic advection in a sim- 
ple kinematic model of a meandering two- 
dimensional jet flow with a Bickley zonal ve- 
locity profile u ~ sech^ y being motivated by 
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the problem of lateral mixing of water masses 
(together with salinity, heat, nutrients, pollu- 
tants, and other passive scalars) in the west- 
ern boundary currents in the ocean. We de- 
rived advection equations in the frame mov- 
ing with the phase velocity of a running wave 
imposed on the Bickley jet, found the station- 
ary points, conditions of their stability, and 
all the possible bifurcations of these equations 
which were shown to be autonomous in the 
co-moving frame. Under a periodic perturba- 
tion of the wave amplitude, the phase plane 
of the chosen model flow has been shown 
to consist of a central eastward jet, periph- 
eral westward currents, and chains of north- 
ern and southern circulations (vortex cores) 
immersed in a chaotic sea which, in turn, 
contains islands of regular motion. Statis- 
tical properties of chaotic transport of ad- 
vected particles have been characterized in 
terms of particle's zonal flights (any event 
between two successive changes of the sign 
of the particle's zonal velocity u). Probabil- 
ity density functions (PDFs) of durations and 
lengths of flights, computed with a number 
of very long chaotic trajectories, have been 
found to be complicated functions with local 
maxima and fragments with exponential and 
power-law decays. 

The aim of this paper is to study a phe- 
nomenological connection between dynami- 
cal, topological, and statistical properties of 
chaotic mixing and transport in the mean- 
dering jet flow considered in Paper I and to 
explain transport properties by phenomenon 
of the so-called dynamical traps. Following 
to Zaslavsky [2], the dynamical trap is a do- 
main in the phase space of a Hamiltonian sys- 



tem where a particle (or, its trajectory) can 
spend arbitrary long finite time, performing 
almost regular motion, despite the fact that 
the full trajectory is chaotic in any appropri- 
ate sense. In fact, it is the definition of a 
quasi-trap. Absolute traps, where particles 
could spend an infinite time, are possible in 
Hamiltonian systems only with a zero mea- 
sure set. The dynamical traps are due to 
a stickiness of trajectories to some singular 
domains in the phase space, largely, to the 
boundaries of resonant islands, saddle trajec- 
tories, and cantori. There are no classifica- 
tion and description of the dynamical traps. 
Zaslavsky described two types of dynamical 
traps in Hamiltonian systems: hierarchical- 
islands traps around chains of resonant is- 
lands f27\ 112] and stochastic-layer traps 
which are stochastic jets inside a stochastic 
sea where trajectories can spend a very long 
time [29| |2l [12]. It is expected that classi- 
fication and description of the most typical 
dynamical traps would help us to construct 
kinetic equations which will be able to de- 
scribe transport properties of chaotic systems 
including anomalous ones j2l [111 [12] . 

This paper is organized as follows. We 
start with advection equations derived in Pa- 
per I in the frame of reference moving with 
the phase velocity of a meander whose am- 
plitude changes in time periodically. We 
compute in Sec. II PDFs for lengths Xf and 
durations Tf of zonal flights for a number 
of chaotic trajectories and show analytically 
that all the flights start and finish only inside 
a strip confined by two curves whose form 
is defined by the condition u = 0. Some 
prominent peaks in statistics of short flights 
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{\xf\ < 2tt) are proved to be caused by sticki- 
ness of trajectories to the boundaries of rota- 
tional resonant islands filled with regular par- 
ticles rotating in the same frame (Sec. III). 
We call this type of dynamical traps 
rotational-island trap (RIT). In Sec. IV we 
study dynamical traps connected with peri- 
odic saddle trajectories, emerged from sad- 
dle points of the unperturbed system (jl]) un- 
der the perturbation ([5]), and prove numeri- 
cally that the saddle traps (STs) contribute 
to the statistics of short flights as well. An- 
other type of islands, ballistic islands (filled 
with regular particles moving from frame to 
frame), is proved to contribute to the statis- 
tics of long flights (|a;/| ^ 27r) in Sec. V. Both 
the ballistic-islands trap (BITs) and RITs be- 
long to the class of hierarchical-island traps 
by the Zaslavsky's classification. 

2 Description of chaotic 
transport in terms of 
flights 

2.1 Basic features of the flow 

We take the following specific stream function 
as a kinematic model of a meandering jet flow 
in the laboratory frame of reference: 



-vl^gtanh 



y' — a cos k{x' — ct') 
Av^l + W sm^k{x' - ct' 



plitude a, the wave-number k, and the phase 
velocity c. The normalized streamfunction in 
the frame moving with the phase velocity is 



\1/ = - tanh 



y — A cos X 
LVI + sin^ 



X 



Cy, (3) 



where x = k{x' — ct') and y = ky' are new 
scaled coordinates, and A = ak, L = Xk, 
and C = c/\I'qA; are the control parame- 
ters. Equations, governing advection of pas- 
sive particles ([1]) in the co-moving frame, are 
the following: 



X 



LVI + A"^ sin^x cosh^ 9 
A sin x(l + A"^ — Ay cos x) 

L{l + A^ sin^ xY^^osh^d'' 
y — A cos X 



(4) 



e 



lVi + A^ 



^2 ' 

sm X 



(2) 

where the width of the jet is A. Meandering 
is provided by a running wave with the am- 



where dot denotes differentiation with re- 
spect to dimensionless time t = "^'^k^t'. 

In Paper I (for more details see |T7]) we 
have found and analyzed all the stationary 
points, their stability, and the bifurcations 
of the equations of motion Being moti- 
vated by the problem of mixing and trans- 
port of water masses and their properties 
in oceanic western boundary currents like 
the Gulf Stream and the Kuroshio, we chose 
the phase portrait shown in Fig. [T^ among 
all the possible flow regimes. Passive par- 
ticles can move along stationary (in the co- 
moving frame) streamlines in a different man- 
ner. They can move to the east in the jet 
(J) and to the west in northern and southern 
(with respect to the jet) peripheral currents 
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(P). There are also particles rotating in the 
northern and southern circulation cells (C) 
in a periodic way. The northern separatrix 
connects the saddle points at a;i"'* = 2?™ and 
j/i"'' = L Arcosh ^JXjLC + A and the south- 
ern one connects the saddle points at xf^ = 
{2n + l)n and yi"^ = -L Arcosh ^l/LC - A, 
where = 0, ±1, .... 

As a perturbation, we took in Paper I the 
simple periodic modulation of the meander's 
amplitude 

A = Aq + e cos{ujt + (p) . (5) 

Under the perturbation, there arise reso- 
nances between the perturbation frequency u 
and the frequencies / of the particle's rota- 
tion in the circulations C. A frequency map 
/(xo,?/o); computed in Paper I (see Fig. 2 in 
that paper), shows the values of / for parti- 
cles with initial positions {xo,yo) in the un- 
perturbed flow. With a given value of the 
perturbation frequency and fixed values of 
the other control parameters, the vortex cores 
in the circulations survive, stochastic layers 
appear along the unperturbed separatrix, and 
the central jet J is a barrier to transport of 
particles across the jet. In Paper I we fixed 
the scaled values of the parameters of the un- 
perturbed flow, the jet's width L = 0.628, 
the meander's amplitude Aq = 0.785 and 
its phase velocity C = 0.1168, that are in 
the range of the realistic values for the Gulf 
Stream [221 |23], and took the initial phase 
to be (j) = 7r/2. The perturbation frequency 
u = 0.2536 chosen in Paper I is close to the 
values of the rotation frequency / of the par- 
ticles circulating in the inner core of the re- 
gions C (see Fig. 2 in Paper I). In Fig. [T]b we 



show the Poincare section (for a large number 
of trajectories) of the meandering jet whose 
amplitude is modulated with the frequency 
w = 0.2536 and the strength e = 0.0785. 
The vortex cores survive under this pertur- 
bation, the stochastic layers appear along the 
unperturbed separatrix, and a central jet J is 
a barrier to transport of particles across the 
jet. 

The equations of motion (jl]) with the per- 
turbation are symmetric under the follow- 
ing transformations: (l)t— j-t, x— >7r + x, 
y — )■ —y and (2) t — )■ —t, x — )■ — x, y y. It 
implies that the meridional transport (north- 
south and south-north) is symmetric but the 
zonal transport (west-east and east- west) is 
symmetric under a time reversal. Due to 
these symmetries motion can be considered 
on the cylinder with < x < 27r and y > 0. 
The part of the phase space with 2™ < x < 
27r(n + 1), n = 0, ±1, . . . , is called a frame. 

It is convenient to characterize chaotic mix- 
ing and transport in terms of zonal flights. A 
zonal flight is a motion of a particle between 
two successive changes of signs of its zonal 
velocity, i. e. the motion between two succes- 
sive events x = u = 0. Particles (and cor- 
responding trajectories) in chaotic jet flows 
can be classified in terms of the lengths of 
flights X/ as follows. The trajectories with 
|x/| < 27r correspond to the particles moving 
in the same frame or in neighbor frames. In 
the global stochastic layer there are particles 
moving chaotically forever in the same frame 
but they are of a zero measure. Among the 
particles with inter-frame motion, there are 
regular and chaotic ballistic ones. Regular 
ballistic trajectories can be defined as those 
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which cannot have two flights with \xf\ > 27i 
in succession. They correspond to particles 
moving in regular regions of the phase space 
persisting under the perturbation, (eastward 
motion in the jet and western motion in 
the peripheral current) and those moving in 
the stochastic layer (trajectories belonging to 
ballistic islands). Typical chaotic trajecto- 
ries have complicated distributions over the 
lengths and durations of flights. 

In the laboratory frame of reference, all 
the fluid particles move to the east together 
with the jet flow and a flight is a motion be- 
tween two successive events when the parti- 
cle's zonal velocity U is equal to the mean- 
der's phase velocity c. U U < c, the corre- 
sponding particle is left behind the meander 
(it is a western flight in the co-moving frame), 
if U > c, it passes the meander (an eastern 
flight in the co-moving frame). Short flights 
with \xf\ < 271 (motion in the same spatial 
frame in the co-moving frame of reference) 
correspond to the motion in the laboratory 
frame when two successive events U = c oc- 
cur on the space interval less than the me- 
ander's spatial period 2iT/k. Ballistic flights 
between the spatial frames in the co-moving 
frame with \xf\ > 27r correspond to the mo- 
tion in the laboratory frame when the parti- 
cles move through more than one meander's 
crest between two successive events U = c. 

2.2 Turning points 

As in Paper I, we will characterize statisti- 
cal properties of chaotic transport by proba- 
bility density functions (PDFs) of lengths of 
flights P{xf) and durations of flights P{Tf) 



for a number of very long chaotic trajec- 
tory. Both regular and chaotic particles may 
change many times the sign of their zonal ve- 
locity X = u. From the condition i = in 
Eq. (jl]), it is easy to find the equations for 
the curves which are loci of turning points 



Y±{x,A) = ±LV 1 + sin^ XX 



Arsech y LC \/ 1 + A"^ sin^ x + A cos x. 



(6) 



We consider the northern curve, i. e., Eq. ([6]) 
with the positive sign. Taking into account 
that the perturbation has the form ([5]), we 
realize that all the northern turning points 
are inside a strip confined by two curves of 
the form with A = Aq^e. Let us analyze 
the derivative over the varying parameter A 



dY 
dA 

ACL"^ sin^ X 



COSX+ 

H I 2 Arsech VD 



2D 



(7) 



where D = LCVl + A^ sin^ X. If the deriva- 
tive at a fixed value of x does not change its 
sign on the interval Aq — e < A < Aq + e, then 
Y varies from Y{x, Aq — e) to Y{x, + e), 
and for each value of y we have a single value 
of the perturbation parameter A. However, 
there may exist such values of x for which 
the equation dY/dA = has a solution on 
the interval mentioned above. In this case 
one may have more than one values of A for 
a single value of y. Thus, the width of the 
strip, containing turning points, is defined by 
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the values of Y at the extremum points and 
at the end points of the interval of the val- 
ues of A. In Fig. [1] c we show the turning 
points of a single chaotic trajectory on the 
cylinder < x < 2tt confined between two 
corresponding curves. 

In the numerical simulation throughout the 
paper we use the Runge-Kutta integration 
scheme of the fourth order with the constant 
time step AtO.0247. To study chaotic trans- 
port we have carried out numerical exper- 
iments with tracers initially placed in the 
stochastic layer. It was found that statisti- 
cal properties of chaotic transport practically 
do not depend on the number of tracers pro- 
vided that the corresponding trajectories are 
sufficiently long {t ~ 10^). The PDFs for the 
lengths Xf and durations Tf of flights for five 
tracers with the computation time t = 5 ■ 10^ 
for each tracer are shown in Figs. |2] a and 
b, respectively, both for the eastward (e) and 
westward (w) motion. Both P{xf) and P{Tf) 
are complicated functions with local extrema 
decaying in a different manner for different 
ranges of xj and Tj. The main aim of our 
study of chaotic transport is to figure out 
the basic peculiarities of the statistics and 
attribute them to specific zones in the phase 
space, namely, to dynamical traps strongly 
infiuencing the transport. 

3 Rotational-islands traps 

It is well known, that in nonlinear Hamilto- 
nian systems a complicated structure of the 
phase space with islands, stochastic layers, 
and chains of islands, immersed in a stochas- 



tic sea, arises under a perturbation due to 
a variety of nonlinear resonances and their 
overlapping [24j. The motion is quasiperi- 
odic and stable in the islands. The bound- 
aries of the islands are absolute barriers to 
transport: particles can not go through them 
neither from inside nor from outside. Invari- 
ant curves of the unperturbed system (see 
Fig. [1] a) are destroyed under the perturba- 
tion ([5]) (see Fig. [T]b). As the perturbation 
strength e increases, a closed invariant curve 
with frequency / is destroyed at some critical 
value of e. If the f /uj is a rational number, 
the corresponding curve is replaces by an is- 
land chain, while the curves with irrational 
frequencies are replaced by cantori (for a re- 
view see There are uncountably many 
cantori forming a complicated hierarchy. Nu- 
merical experiments with a variety of Hamil- 
tonian systems with different number of de- 
grees of freedom provide an evidence for the 
presence of strong partial barriers to trans- 
port around the island's boundaries (for re- 
view, see [12]) which manifest themselves on 
Poincare sections as domains with increased 
density of points. 

In Paper I we have found that with chosen 
values of the control parameters there exist 
in each frame a vortex core (which is an is- 
land of the primary resonance u = f) im- 
mersed into a stochastic sea, where there are 
six islands of a secondary resonance emerged 
from three islands of the primary resonance 
3/ = 2lj (see Fig. 3 in Paper I). Chains of 
smaller size islands are present around the 
vortex core and the secondary-resonance is- 
lands. Particles belonging to all of these is- 
lands (including the vortex core) rotate in the 
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same frame performing short flights with the 
lengths \xf\ < 271. So we will call them ro- 
tational islands and distinguish from the so- 
called ballistic islands to be considered below. 

Stickiness of particles to boundaries of the 
rotational islands has been demonstrated in 
Paper I. It means that real fluid particles 
can be trapped for a long time in a singu- 
lar zone nearby the borders of the rotational 
islands which we will call rotational-islands 
traps (RITs). To illustrate the effect of the 
RITs we demonstrate in Figs. [2] and H] the 
Poincare sections of a chaotic trajectory in 
the frame < a; < 27r sticking to the vortex 
core and to the secondary-resonance islands, 
respectively. The contour of the vortex core 
is shown in Fig. H] by the thick line. The 
small points are tracks of the particle's po- 
sition at the moments of time in = 27rn/u 
(where n = 1,2,...) and the thin curves 
are fragments of the corresponding trajec- 
tory on the phase plane. Increased density 
of points indicates the presence of dynamical 
traps near the boundaries of the rotational 
islands. Contribution of the vortex-core RIT 
(Fig. [3]) to chaotic transport is expected to 
be much more significant than the one of the 
RITs of the other islands (Fig. H]). 

It is reasonable to suppose that RITs con- 
tribute to the statistics of short flights. By 
short flights we mean the flights with the 
length shorter than 271. In Fig. |5]we show the 
part of the full PDF P{Tf) (Fig. [2]b) for the 
eastward (e) and westward (w) short flights 
separately. There are a comparatively small 
number of the eastward flights with Tf < 11. 
Let us note the prominent peak of the cor- 
responding PDF at ^ 11 followed by an 



exponential decay. As to the westward short 
flights, there are two small local peaks around 
Tf ~ 17 and 21. 

To estimate the contribution of the vortex- 
core RIT to the statistics of short flights, we 
compute and compare the statistics of the du- 
rations of flights Tf for two trajectories: a 
regular quasiperiodic one with the initial po- 
sition close to the inner border of the vortex 
core (Fig. [6] a) and a chaotic one with the ini- 
tial position close to the vortex-core border 
from the outside (Fig. E]b). Each full rota- 
tion of a particle in a frame consists of two 
flights, eastward and westward, with different 
values of Tf because of the zonal asymmetry 
of the flow. The statistics for the chaotic tra- 
jectory, sticking to the vortex core (Fig. E]b), 
may be considered clS db distribution of the du- 
rations of flights in the vortex-core RIT. The 
minimal flight duration in this RIT is ~ 11 
(the flights with smaller values of Tf are rare 
and they occur outside the trap). Positions 
of the local maxima of the PDF for the stick- 
ing trajectory in Fig. [6] b correlate approxi- 
mately with the corresponding local maxima 
of the PDF for the regular trajectory inside 
the core in Fig. [Hla. The similar correlations 
have been found (but not shown here) be- 
tween the local maxima of the PDFs for the 
lengths of flights P{xf ) for the interior reg- 
ular and sticking chaotic trajectories. These 
correlations and positions of the peaks prove 
numerically that short flights with \xf\ < 27r 
and 11 < T/ < 21 may be caused by the 
effect of vortex-core RIT. We conclude from 
Fig. [6]b that the vortex-core RIT contributes 
to the statistics of the short flights in the 
range 11 < < 20 for the eastward flights 
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with the prominent peak at Tj ^ 11 and in 
the range 15 < < 21 for the westward 
flights with small peaks at Tf ~ 17 and 21. 

The effect of the RIT of the secondary- 
resonance islands is illustrated in Fig. HI 
To find the characteristic times of this RIT 
we compute two trajectories: a regular 
quasiperiodic one with the initial position in- 
side one of these islands and a chaotic one 
with the initial position close to the outer 
border of the island. The respective PDFs 
P{Tf), shown in Figs. [7]a and b, demonstrate 
strong correlations between the correspond- 
ing peaks at Tf ~ 12, 23, and 27. Computed 
(but not shown here) PDFs P{xf) for these 
trajectories confirm the effect of the islands 
RIT on the statistics of short flights. 

4 Saddle traps 

As a result of the periodic perturbation (l5l), 
the saddle points of the unperturbed system 
dH) at xi"^ = 27rn, yi"^ = L Arcosh ^1/LC + 

A and at xj'^ = {2n + l)7r, yi'^ = 
-L Arcosh ^/l/LC - A (n = 0,±1,...) be- 
come periodic saddle trajectories. These hy- 
perbolic trajectories have their own stable 
and unstable manifolds and play a role of 
specific dynamical traps which we call saddle 
traps (ST). In this section we demonstrate 
that the STs influence strongly on chaotic 
mixing and transport of passive particles and 
contribute, mainly, in the short-time statis- 
tics of flights. 

Tracers with initial positions close to a 
stable manifold of a saddle trajectory are 
trapped for a while performing a large num- 



ber of revolutions along it. To illustrate the 
effect of the STs we show in Fig. [H] a and 
b fragments of two chaotic trajectories stick- 
ing to the saddle trajectory and performing 
about 20 full revolutions before escaping to 
the east (Fig. [H]a) and to the west (Fig. |H]b). 
We have managed to detect and locate the 
corresponding periodic unstable saddle tra- 
jectory which is situated in Figs. [8] a and 
b in the domain where a few fragments of 
the chaotic trajectory imposed on each other. 
Because of the flow asymmetry, the dura- 
tion of eastern flights of a particle along the 
saddle trajectory Te ~ 11.9 is shorter than 
the duration of western flights — 12.9. 
The black points are the tracks of the par- 
ticle's positions on the flow plane at the mo- 
ments of time tn = 2im/uj ~ 24.8 n (where 
n = 1, 2, . . . ). They belong to smooth curves 
which are fragments of the stable and unsta- 
ble manifolds of the saddle trajectory at the 
chosen initial phase = tx/2. 

To estimate the contribution of the STs to 
the statistics of short flights shown in Fig. [5l 
we compute and plot in Fig.[8]c the number of 
the eastward (A^e) and westward (A^w) short 
flights with a given duration Tf for those two 
chaotic trajectories sticking to the saddle tra- 
jectory arising from the saddle point with the 
position Xs = 0,ys — 2.02878. Each full ro- 
tation of the particles consists of an eastward 
flight with the duration Tg ~ 11.9 and an 
westward flight with the duration T„ ~ 12.9. 
The flights with Te ~ 11.9 contribute to 
the main peak in Fig. [5] and the flights with 
T„ ~ 12.9 to "the wesward" plateau in that 
figure. 

The mechanism of operation of the STs 



9 



can be described as follows. Each saddle 
trajectory 7(t) possesses time-dependent sta- 
ble Wsijlt)) and unstable Wuijit)) material 
manifolds composed of a continuous sets of 
points through which pass at time t trajec- 
tories of fluid particles that are asymptotic 
to 7(t) as t — > oo and t — )■ — oo, respectively. 
Under a periodic perturbation, the stable and 
unstable manifolds oscillate with the period 
of the perturbation. It was firstly proved by 
Poincare that Ws and may intersect each 
other transversally at an infinite number of 
homoclinic points through which pass doubly 
asymptotic trajectories. To give an image of 
a fragment of the stable manifold of the peri- 
odic saddle trajectory, we distribute homoge- 
neously 2.5 ■ 10^ particles in the rectangular 
[-0.4 < a; < 0.45; 2 <y < 2.1] and compute 
the time the particles need to escape the rect- 
angular. The color in Fig. |9] modulates the 
time T when particles with given initial posi- 
tions (xq, Vq) reach the western line at x = — 1 
or the eastern line at x = 1. The particles 
with initial positions marked by the black and 
white colors move close to the stable manifold 
of the saddle trajectory and spend a maximal 
time near it before escaping. The black and 
white diagonal curve in Fig. is an image of 
a fragment of the corresponding stable man- 
ifold. The particles with initial positions to 
the north from the curve escape to the west 
along the unstable manifold of the saddle tra- 
jectory whereas those with initial positions to 
the south from the curve escape to the east 
along its another unstable manifold. 

We have found that particles quit the ST 
along the unstable manifolds in accordance 
with specific laws. We distribute a large 



number of particles along the segment with 
Xo = and i/q = [2.02; 2.06], crossing the 
stable manifold Wg, and compute the time 
T particles with given initial latitude po- 
sitions Ho need to quit the ST. More pre- 
cisely, T(?/o) is a time moment when a par- 
ticle with the initial position i/q reaches the 
lines with x = — lorx = l. The "experi- 
mental" points in Fig. [10] a fit the law Te = 
(-85.81±0.04)-(31.216±0.007) ln(yo«-l/o) 
for the particles which quit the trap moving 
to the east and the law T„ = (— 60.61±0.03) — 
(28.933 ± 0.006) ln(?/o - Vos) for those parti- 
cles which move to the west when quitting the 
trap, where yos = 2.0405755472 is a crossing 
point of Ws with the segment of initial posi- 
tions. 

The ST attracts particles and force them 
to rotate in its zone of influence performing 
short flights, the number of which n depends 
on particle's initial positions yo. The n(?/o) 
is a steplike function (see Fig. [TU] b) with 
the lengths of the steps decreasing in a ge- 
ometric progression in the direction to the 
singular point, Ij = loq~\ where Ij is the 
length of the j-th step and q ^ 2.27 for the 
western exits and q ~ 2.20 for the eastern 
ones. The seeming deviation from this law 
in the range yo = [2.045; 2.046] (see a small 
western segment between two larger ones in 
Fig. [TU] b) is explained by crossing the ini- 
tial line yo = [2.02; 2.06] by the curve of zero 
zonal velocity u. To have the correct law for 
the western exits, it is necessary to add the 
two segments of that cut step. The asymme- 
try of the functions T(?/o) and n(yo) is caused 
by the asymmetry of the flow. 
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5 Ballistic-islands traps 

Besides the rotational islands with parti- 
cles moving around the corresponding ellip- 
tic points in the same frame, we have found 
in Paper I ballistic islands situated both in 
the stochastic layer and in the peripheral 
currents. Regular ballistic modes [30j corre- 
spond to stable quasiperiodic inter-frame mo- 
tion of particles. Only the ballistic islands in 
the stochastic layer are important for chaotic 
transport. Mapping positions of the regular 
ballistic trajectories at the moments of time 
tn = 2im/u} {n = 1,2,...) onto the first 
frame, we obtain chains of ballistic islands 
both in the northern and southern stochastic 
layers, i. e., between the borders of the north- 
ern (southern) peripheral currents and of the 
corresponding vortex cores. A chain with 
three large ballistic islands is situated in those 
stochastic layers. The particles, belonging to 
these islands, move to the west, and their 
mean zonal velocity can be easily calculated 
to be {uf) = -27r/3T = -u/3 ~ -0.0845. 
There are also chains of smaller-size ballistic 
islands along the very border with the periph- 
eral currents. 

We have demonstrated in Paper I a sticki- 
ness of chaotic trajectories to the borders of 
those three large ballistic islands (see Figs. 6 
and 7 in Paper I). The Poincare section 
with fragments of two chaotic trajectories 
in the northern stochastic layer is shown in 
Fig. [11] a. One particle performs a long flight 
sticking to the very border with the regular 
westward current, and another one moves to 
the west sticking to the very boundaries of 
three large ballistic islands. A magnification 



of a fragment of the border and tracks of a 
sticking trajectory around a smaller-size bal- 
listic island are demonstrated in Fig. [11] b. 
Fig. [11] c demonstrates the effective size of the 
trap of the large ballistic islands with tracks 
of a sticking trajectory around them. 

It is reasonable to suppose that the 
ballistic-islands traps (BIT) contribute, 
largely, to the statistics of long flights 
with \xf\ ^ 2tt. All the ballistic particles, 
moving both to the west and to the east, 
can finish a flight and make a turn only 
in the strip shown in Fig. [1] c. The loci 
of the corresponding turning points have 
a complicated fractal-like structure. We 
consider further only long westward flights, 
taking place in the northern stochastic layer, 
because it is much wider than the stochastic 
layer between the regular central jet and 
the southern parts of the vortex cores where 
eastward flights take place. 

To distinguish between contributions of 
the traps of different ballistic islands (and, 
maybe, other zones in the phase space) to the 
statistics of long flights, we compute for five 
long chaotic trajectories (up tot = 5-10^) the 
distribution of a number of westward flights 
with Tf > 10^ over the mean zonal veloci- 
ties (uf) = Xf/Tf of the particles perform- 
ing such flights. The distribution in Fig. [12] 
has a prominent peak centered at the mean 
zonal velocity (uf) ^ —0.0845 which corre- 
sponds to a large number of long flights of 
those particles (and their trajectories) which 
stick to the very boundaries of the large bal- 
listic islands (see Fig. [11] a) moving with the 
mean velocity {uf) ~ -0.0845. The flat left 
wing of the distribution N{{uf)) corresponds 
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to the traps of smaller-size ballistic islands 
nearby the border with the peripheral cur- 
rent. There are different families of these is- 
lands (see one of them in Fig.[TT]b) with their 
own values of the mean zonal velocity which 
are in the range -0.092 < {uf) < -0.0845. 
Stickiness to the boundaries of the border is- 
lands is weaker because they are smaller than 
the large islands and their contribution to 
the statistics of long flights is comparatively 
small. 

The right wing of the distribution N{{uf)) 
with -0.084 < {uf) < -0.075 deserves fur- 
ther investigation. The value (uf) ~ —0.075 
is a minimal value of the zonal velocity for 
long westward flights possible in the north- 
ern stochastic layer. Increasing the mini- 
mal duration of a flight from Tf = 10^ to 
Tf = {2 5) ■ 10^, we have found splitting of 
the broad distribution with —0.084 ^ (m/) ^ 
—0.08 into a number of small distinct peaks. 
Comparing trajectories with the values of 
(uf) corresponding to these peaks, we have 
found that all they move around the large 
ballistic islands. The particles with smaller 
values of (uf) used to penetrate further to 
the south from the islands more frequently 
than those with larger values of (uf) which 
prefer to spend more time in the northern 
part of the dynamical trap connected with 
those islands. Thus, we attribute the right 
wing of the distribution N{{uf)) to an effect 
of the trap situated around the large ballistic 
islands. 

To estimate the contribution of different 
BITs to the statistics of long westward flights 
in Fig. [2] we have computed the PDFs P{xf) 
and P{Tf) for particles performing westward 



flights with Xf > 100 and Tf > 1000 and 
with the mean zonal velocity (uf) to be cho- 
sen in three different ranges shown in Fig. [121 
-0.092 < (uf) < -0.085 (particles stick- 
ing to the border islands) —0.085 < {uf) < 
—0.084 (particles sticking to the very bound- 
ary of three large islands), and —0.084 < 
(■^z) ~ —0.075 (the trap of the three large 
islands). All the PDFs P{xf) decay ex- 
ponentially but with different values of the 
exponents equal to z/ ^ —0.005 and u ^ 
-0.0018-^ — 0.0014 for the traps of border and 
the large ballistic islands, respectively. The 
tail of the PDF P{xf) for westward flights, 
shown in Fig. [2l decays exponentially with 
u —0.0014. Thus, the contribution of the 
large island's BIT to the statistics of long 
westward flights is dominant. As to tem- 
poral PDFs P{Tf) for westward long flights, 
they are neither exponential nor power-law 
like with strong oscillations at the very tails. 
The slope for the border BITs is again smaller 
than for the large ballistic islands trap. 

6 Conclusion 

A meandering jet is a fundamental structure 
in oceanic and atmospheric flows. We de- 
scribed statistical properties of chaotic mix- 
ing and transport of passive particles in a 
kinematic model of a meandering jet flow in 
terms of dynamical traps in the phase (phys- 
ical) space. The boundaries of rotational is- 
lands (including the vortex cores) in circula- 
tion zones are dynamical traps (RITs) con- 
tributing, mainly, to the statistics of short 
flights with \xf \ < 27r. Characteristic times 
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and spatial scales of the RITs have been 
shown to correlate with the PDFs for the 
lengths Xf and durations Tf of short flights. 
The stable manifolds of periodic saddle tra- 
jectories play a role of saddle traps (STs) with 
the specific values of the lengths and dura- 
tions of short flights of the particles sticking 
to the saddle trajectories. The boundaries of 
ballistic islands in the stochastic layers (in- 
cluding those situated along the border with 
the peripheral current) are dynamical traps 
(BITs) contributing, mainly, to the statistics 
of very long flights with 3> 2tt. 

Dynamical traps are robust structures in 
the phase space of dynamical systems in the 
sense that they present at practically all val- 
ues of the corresponding control parameters. 
We never know exact values of the param- 
eters in real flows, especially, in geophysical 
ones. We don't know exactly the structure 
of the corresponding phase space, however, 
we know that typical features, like islands of 
regular motion, vortices, and jets, exist in 
real flows (see their images in some labora- 
tory flows in Ref. [31j). In this paper we 
chose specific values of the control parame- 
ters for which specific PDFs have been com- 
puted and explained by the effect of those dy- 
namical traps that exist under the chosen pa- 
rameters. We have carried out computer ex- 
periments with different values of the control 
parameters and found that the phase space 
structure has been changed, of course, with 
changing the values of the parameters, but 
the corresponding RITs, STs, and BITs with 
specific temporal and spatial characteristics 
have been found to contribute to the corre- 
sponding statistics. 



After finishing our work, we were ac- 
quainted with Ref. [32j where meridional 
chaotic transport, associated with a similar 
kinematic model of a meandering jet, has 
been studied by the method of lobe dynam- 
ics [33]. In difference from our study of zonal 
chaotic transport, a geometric structure of 
cross jet transport has been considered in 
Ref. [32] where values of the control parame- 
ters have been chosen to be sufficiently large 
to break up the central jet as a barrier to 
transport of particles across the jet. The 
mechanisms for particles to cross the jet have 
been described in terms of lobe dynamics and 
the mean time to cross the jet for particles 
entering the jet and the mean residence time 
for particles in the jet have been estimated 
in Ref. [32]. We have studied a more real- 
istic situation (at least, for surface oceanic 
jet currents) when the jet is an absolute bar- 
rier to cross jet transport and we explained 
statistical properties of transport in terms of 
dynamical traps of saddle trajectories, rota- 
tional and ballistic islands. The method of 
lobe dynamics is hardly applicable for study 
zonal chaotic transport since it is practically 
impossible to trace out lobe evolution for a 
large number of frames. 
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Figure 1: (a) Stationary streamfunction of 
a meandering jet in the co-moving frame 
([2]). The flow is divided into three different 
regimes: circulations (C), jet (J), and pe- 
ripheral currents (P). (b) Poincare section 
of the perturbed meandering jet in the co- 
moving frame. The parameters of the steady 
flow are: the jet's width L = 0.628, the me- 
ander's amplitude Aq = 0.785 and its phase 
velocity C = 0.1168. The perturbation am- 
plitude and frequency are: e = 0.0785 and 
u = 0.2536. (c) Turning points of a single 
chaotic trajectory on the cylinder < x < 271 
are in a strip confined by two curves with 
A = Ao±e. 
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Figure 2: Probability density functions of 
(a) lengths Xf and (b) durations Tf of the 
westward (w) and eastward (e) flights. The 
PDFs P^{Tf) and Pc{Tf) are normahzed to 
the number of westward (4.23 ■ 10'') and east- 
ward (4 • 10^) flights, respectively. Statistics 
for five tracers with the computation time 
t — b ■ 10^ for each one. 
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Figure 3: The vortex-core trap. Poincare 
section of a chaotic trajectory in the frame 
< X < 27r with a fragment of a trajectory. 
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Figure 4: The secondary resonance islands 
trap. A fragment of a chaotic trajectory 
sticking to the islands is shown. 
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Figure 5: The PDFs for the eastward (e) and 
westward (w) flights with the length shorter 
than 27r. The PDFs P^{Tf) and Pe{Tf) 
are normahzed to the number of westward 
(4.19 ■ 10^) and eastward (3.7 ■ 10^) flights, 
respectively. Statistics for five tracers with 
the computation time t — 10^ for each one. 
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Figure 6: The vortex-core trap PDFs of du- 
rations Tf of the eastward (e) and westward 
(w) flights, (a) Regular quasiperiodic trajec- 
tory with the duration t = 2 ■ 10^ inside the 
vortex core close to its boundary. Both the 
PDFs are normalized to the number 8 ■ 10^ of 
corresponding flights, (b) Chaotic trajectory 
with the duration t — 2 ■ 10^ sticking to the 
boundary of the vortex core from the outside. 
Both the PDFs are normalized to the number 
4 • 10^ of corresponding flights. 
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Figure 7: The secondary- resonance islands 
trap. The PDFs of durations Tf of the east- 
ward (c) and westward (w) flights, (a) Regu- 
lar quasiperiodic trajectory inside the islands 
with the duration i = 5 • 10^ Both the PDFs 
are normalized to the number 1.5T0^ of corre- 
sponding flights, (b) Chaotic trajectory stick- 
ing to the island's boundary from the outside 
with the duration t = 5 ■ 10^. P^{Tf) and 
Pe(Tf) are normalized to the number of west- 
ward (1.1 • 10^) and eastward (9 • 10^) flights, 
respectively. 
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Figure 8: The saddle trap. Fragments of two 

chaotic trajectories sticking to the periodic 
saddle trajectory one of which escapes to the 
east (a) and another one to the west (b). (c) 
The number of the eastward (iVe) and west- 
ward (A^w) short flights with duration Tf for 
those two trajectories. Statistics with two 
trajectories with the duration t = 10'^ and 
the total number of western = 55 and 
eastern A^e = 51 flights. 



Figure 9: The saddle-trap map. Color modu- 
lates the time T which 2.5-10^ particles with 
given initial positions {xq, Uq) need to reach 
the lines ata; = — lora; = l escaping to 
the west (w) and to the east (e), respectively. 
The black and white diagonal curve is an im- 
age of a fragment of the stable manifold of the 
saddle trajectory. The cross is a position of a 
particle on that trajectory at the initial time 
moment. The integration time is i = 500. 
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Figure 10: (a) Time T a particle with an ini- 
tial latitude position i/q needs to quit the sad- 
dle trap, (b) The number of short flights n 
such a particle performs before quitting the 
saddle trap. The ranges of t/o from which par- 
ticles quit the trap moving to the west and 
east are denoted by "w" and "e", respectively. 
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Figure 11: (a) Poincare section of the north- 
ern stochastic layer where stickiness to the 
very border with the regular westward cur- 
rent and to three large ballistic islands are 
shown. Increased density of points along the 
border with the peripheral current is caused 
by the traps of the border baUistic islands one 
of which is shown in (b) . (c) The trap of the 
large ballistic islands. 
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Figure 12: The distribution of a number of 
long westward flights with Tf > 10'^ over their 
mean zonal velocities (uf). The sharp peak 
corresponds to the trap connected with the 
very boundaries of the large baUistic islands, 
the left wing — to a number of traps of fam- 
ilies of the border ballistic islands, and the 
right wing — to the trap situated around the 
large baUistic islands. Statistics for five trac- 
ers with the total number of long westward 
flights A^j = 5 ■ 10^ and the computation time 
i = 5 • 10^ for each tracer. 
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